New gapped quantum phases for S = 1 spin chain with D2h symmetry 
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We study different quantum pliases in integer spin systems with on-site D2h ~ D2 ® Z2 and 
translation symmetry. We find four distinct non-trivial phases in S = 1 spin chains despite they all 
have the same symmetry. AH the four phases have gapped bulk excitations, doubly-degenerate end 
states and the doubly-degenerate entanglement spectrum. These non-trivial phases are examples of 
symmetry protected topological (SPT) phases introduced by Gu and Wen. One of the SPT phase 
correspond to the Haldane phase and the other three are new. These four SPT phases can be 
distinguished experimentally by their different response of the end states to weak external magnetic 
fields. According to Chen-Gu-Wen classification, the D2h symmetric spin chain can have totally 64 
SPT phases that do not break the symmetry. Here we constructed seven nontrivial phases from the 
seven classes of nontrivial projective representations of D2h group. Four of these are found in S = 1 
spin chains and studied in this paper, and the other three may be realized in 5* = 1 spin ladders or 
S — 2 models. 
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I. INTRODUCTION 

Topological order was introduced to distinguish dif- 
ferent phases which can not be separated by symmetry 
breaking orders.^ Using a definition of phase and phase 
transition based on local unitary transformations, Ref. |3 
shows that what topological order really describes is ac- 
tually the pattern of long range entanglements in gapped 
quantum systems. 

The Haldane phase'^ for S" = 1 spin chains was regarded 
as a simple example of topological order. For a long 
time, the existence of string order (or hidden Z2 Z2 
symmetry breaking), nearly degenerate end states and 
gapped excitations were considered as the hallmark of the 
Haldane phase.'* However, it was shown that even after 
we break the spin rotation symmetry which destroy the 
string order and gapped end states, the Haldane phase 
can still exist (i.e. is still distinct from the trivial phase). 
It was also shown that the Haldane phase has no long 
range entanglements.^^ In fact, all ID gapped ground 
state has no long range entanglements.^ Thus there are 
no intrinsic topologically ordered states in gapped ID 
systems^ (including the Haldane phase). 

But when Hamiltonians have some symmetries, even 
short range entangled states with the same symmetry can 
belong to different phases.^''^ Such phases are called 'sym- 
metry protected topological (SPT) phases' by Gu and 
Weni^ We would like to remark that, due to their triv- 
ial intrinsic topological order and short range entangle- 
ments, 'symmetry protected topological phases' should 
be more properly referred as 'symmetry protected short 
range entangled phases'. 

Thus the Haldane phase is not an intrinsically topo- 
logically ordered phase, but actually an example of SPT 
phase protected by translation and SO (3) spin rotation 
symmetries. This result is supported by a recent real- 
ization that the existence of the Haldane phase requires 
symmetry (such as parity, time reversal, or spin rota- 



tional symmetry) j^i^ In other words, if the necessary sym- 
metries are absent, the Haldane phase can continuously 
connect to the trivial phase without any phase transi- 
tion. We would like to mention that, strictly speaking, 
the topological insulatorsiSrJ^ are not intrinsically topo- 
logicaly ordered phases either. They are other examples 
of SPT phases protected by time reversal symmetry. 

Later, Haldane and his collaborators discovered that 
the entanglement spectrum of the Haldane phase (and 
other nontrivial phases) is doubly degenerate.^ The en- 
tanglement spectrum degeneracy is symmetry-protected. 
A systematic classification of all SPT orders for all ID 
gapped systems was obtained in Ref. ji, using projective 
representations. All ID gapped phases are either sym- 
metry breaking phases or SPT phases. 
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FIG. 1: (Color online) The phase diagram of model ([T}. The 
transition between the Neel phase and the SPT phases are 
second order, and the transition between To and is first 
order. 

In this paper, we will study the SPT phases of a spin- 
chain protected by translational symmetry and on-site 
-D2/1 = D2 <E) Z2 symmetry, where the Z2 in D2 Z2 is 
generated by the time reversal symmetry T. Here we as- 
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sume that the physical spin forms a hnear representation 
of D2 and — 1. First, we study a simple spin-1 model 
with those symmetries 

H = ^^cos9Sa:,iSx,i+i+sm9[cos4){Sy^iSy^i+i (1) 

i 

-^Sz,iSz,i+l) + sin (j){Sxz,iSxz,i+l + Sxy.iSxy.i+l)] ■ 

where Smn = SmSn + SnSm {m, n = x,y, z). As shown 
in Fig. [U this model has three phases, the Neel phase, 
the To phase and the Tx phase. The Neel phase breaks 
the D2h symmetry, and the other two phases do not 
break any symmetry. We can use sublattice spin mag- 
netization as an order parameter to distinguish the Neel 
phase. However, the remaining two phases cannot be dis- 
tinguished through local order parameter such as sublat- 
tice spin magnetization. Further, in the both phases, the 
entanglement spectrum- is doubly degenerate, so both 
of them are nontrivial. However, the entanglement spec- 
trum is not a good order parameter to separate them. 
Therefore, we need a new tool to distinguish these two 
non-trivial phases that cannot be described by symmetry 
breaking and the entanglement spectrum. It turns out 
that the new tool is the projective representation of the 
symmetry group. The two phases can be distinguished 
since their doubly degenerate end states form different 
projective representations of D2h- The physical conse- 
quence is that they response differently to weak external 
magnetic field. The doubly degenerate end states can 
be viewed as an effective spin- 1/2 spin with asymmet- 
ric g-factors: Qx and gy and Qz describing the coupling 
of the end spin to external magnetic field in x, y and z 
directions. We find that gx,gy,gz 7^ in the Tp-phase 
and gx ^ Q,gy = gz = ^ in the Ta;-phase. We would 
like to stress that such a property is robust against any 
perturbations that do not break the D2h symmetry (the 
perturbation may even break the translation symmetry). 

The D2h symmetric spin-chain can have very rich 
quantum phases. It is shown that it can have 64 different 
gapped phases that do not break the Z?2/i and the trans- 
lation symmetry. In fact, it is the projective representa- 
tion theory that allow one to find all the non-trivial SPT 
phases beyond the symmetry breaking description. In 
this paper, we will not study all of them. We will only use 
8 classes of projective representations of the D2h group 
to construct 8 gapped no-symmetry-breaking phases, one 
is trivial and the other 7 are nontrivial SPT phases. We 
find that four of the 7 SPT phases (labeled as Tq, X,, Ty 
and Tz) can be realized in 5 = 1 spin chain. Here Tq is 
the usual Haldane phase (because it includes the Heisen- 
berg point), and Tx,Ty,Tz are new SPT phases. The 
states in different SPT phases cannot be smoothly con- 
nected to each other without explicitly breaking the D2h 
symmetry in the Hamiltonian. The remaining three SPT 
phases cannot be realized for 5=1 chains, but they may 
exist in 5 = 1 spin ladders or 5* = 2 spin chains. 

The four SPT phases are experimentally distinguish- 
able due to the different behaviors of their end states. In 



the To phase, the end states can be considered as spin- 
1/2 free spins. So weak magnetic field couples to the 
end spins and lifts the ground state degeneracy at lin- 
ear order. However, in Tx phase, the end states can no 
longer be considered as normal spin- 1/2 spins because 
they behave differently under time reversal. As men- 
tioned above, the g-factors gy,gz ~ 0, which means that 
By and Bz can not split the degenerate ground states 
in Tx phase at linear order. Similarly, the end states of 
Ty (or Tz) only respond to By (or Bz). According to 
these properties, we propose an experimental scenario to 
distinguish these four phases. 

This paper is organized as the following. In section 
II, we introduce the four SPT phases for the S = 1 spin 
chain models. In section HI, we focus on the interaction 
of the end states to weak external magnetic fields and 
propose an experimental method to distinguish different 
SPT phases. In section IV we briefly summarize the re- 
lationship between the SPT phases and the classes of 
projective representations, and leave detailed derivations 
to the appendix. Section V is devoted to conclusions and 
discussions. 



II. THE MODEL AND SPT PHASES 

D2h group has eight group elements, D2h — 
{E,Rx,Ry,Rz,T,RxT,RyT,RzT}, which is a direct 
product of the 180° spin rotation group D2 — {E, Rx — 

g-i7r5^ , Ry — e^'i-^Sy ^ _ g-j7r5j j. g^j^j iime reversal 

symmetry group Z2 = {i?,T}. Note that T inverts the 
spin [Sx, Sy, Sz) — > {^Sx, ~Sy, —Sz) and is anti-unitary. 
D2h has eight 1-d linear representations (as shown in 
Tab. lVIl in appendix lBj) . Since T is anti- unitary, the bases 
|(/)) and i\<j)) have different time reversal parity. This sub- 
tle property yields more than one SPT phases. 

The most general Hamiltonian for an 5 = 1 spin chain 
with D2h symmetry and with only nearest neighbor in- 
teraction is given by 

HD2h — ''^^[^i^'xA^'x.j + '^^{Sx.iSy J + Syj^Sxj) 

i 

+a3Sy^Syj + a4{S'^^S'^ J + S'^^^S'^j) 
+a5Sz,iSzj + o,6{Sy^iSzj + Sz,iSyj) 

~^C2Sxz,iSxz,j H~ diSz^iSz.j + d2Sxy,iSxyJ 

+eiSl^ + e2Sl^ + e3Sl^]. (2) 

where j = i + I, Smn = SmSn + SnSm im,n = x,y,z) 
and fli, a2, ei, 62, 63 are constants. We are interested 
in the parameter regions within which the excitations are 
gapped and the ground states respect the D2h symmetry. 

In general, for ID systems with translation symmetry 
and on-site symmetry group G, a gapped ground state 
that does not break any symmetry can be approximately 
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written as a matrix product state (MPS) 



we obtain 



J2 Tr(A™i...A™")|mi...mjv), 



(3) 



{mi,...,mjv} 



which varies in the following way under the symmetry 
group 



(4) 



where g £ G is a group element, and a{g)/M{g) is 
its linear /projective representation matrix. Thus it is 
concluded that the SPT phases are classified by (a;,a), 
where w is the element of the second cohomology group 
H^{G,C) (which describe different classes of projective 
representations of the symmetry group G).^ 

So in our case, the ground state can be generally writ- 
ten in forms of MPS shown in Eq. ([3]) . The requirement 
\ip) being invariant under D2h is equivalent to the condi- 
tion in Eq. (|4|). The main task of this paper is to try to 
find different kinds of states that satisfy this condition. 
In this paper, we only consider the case a{g) = 1. The 
full classification (with a different approach from this pa- 
per) is given in Ref. [IB- 
Let us first consider the on-site terms in ©. When 
|ei|, 1 62 1 or 1 63 1 is large, the ground state of Hu^^ is sim- 
ple. For instance, when 63 — > —00, the ground state is a 
long-range ordered state which breaks the D2h symme- 
try; when 63 — cx), the ground state is a product state 
= Yii (8)10)4, which is trivial Since we are interested 
in the nontrivial SPT phases, we will set ei = 62 = 63 = 
in the following discussion. 



A. Exactly solvable models 

The Affleck Kennedy Lieb Tasaki (AKLT) model^^ is 
an exactly solvable model with 50(3) symmetry that 
falls in the Haldane phase. The AKLT model contains 
all the physical properties of the Haldane phase and all 
other states in this phase can smoothly deform into the 
AKLT state. Since the ground state wave function of this 
exactly solvable model is a simple matrix product (sMP) 
statei^ and is known in advance, so studying this model 
is relatively easy and helps to understand the physics of 
the Haldane phase. 

In this section, we will introduce four classes of exactly 
solvable models that have D2h symmetry. Analogous to 
the AKLT state, the ground states of these exactly solv- 
able models are nontrivial sMP states satisfying Eq. (U). 
We will show that different classes of sMP states can not 
be smoothly connected, which indicates that each class 
corresponds to a phase. 

The first example is a direct generalization of the 
AKLT model. The ground state of AKLT state is rep- 
resented by = (7j;, = ay, — a^, which has 
SO{3) symmetry. When generalized to D2h symmetry. 



aax, A^ = bay, A^ = ca^ 



(5) 



where a, 6, c are nonzero real numbers (the same below). 
When a = 6 = c = 1, the above state reduces to the 
AKLT state. For this reason, we say that this model also 
belongs to the Haldane phase. We label this phase as 
Tq. Similar to the AKLT model, the parent Hamiltonian 
of above state is composed by projectors (for details see 
appendix El and [B]) 



^0 = E 



+(i + a'^b'^j)Sz^iSz,j + - b'^c^l)Syz,iSyz,j 

+ (^ - a^C^'l)Sxz,iSxz,j + " a^h'^l)Sxy,iSxy,j 



ho. 



(6) 



and 



2 

2 J 



where 7 = 2(,4+fc4+,4) 

i 



At open boundary condition, the Hamiltonian ^ has 
exactly four-fold degenerate ground states independent 
of the chain length. The above exactly solvable model 
is frustration free, that is, the expectation value of the 
Hamiltonian is minimized locally in the ground states. 
The excitations are gapped and all correlation functions 
of local operators are short ranged. Furthermore, if a, h, c 
are normalized + b'^ + = 1, then it is easily checked 
that 



here A = /, indicating that the entanglement spectrum 
of the ground states is doubly degenerate. This informs 
that state ^ is nontrivial. Actually, the models at the 
vicinity of ([6]) (the phase To) have very similar properties 
unless gap closing (second order phase transition) or level 
crossing (first order phase transition) happens. 

Now we consider another example of sMP state. 



A"" = 



A^ = bay, A" = cc 



(7) 



Above sMP state is also invariant under D2h group. As 
will be shown later, it can not be continuously connected 
to Eq. without breaking the D2h symmetry. This 
means that it belongs to another phase which we label 
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as Tx phase. The parent Hamihonian of ([7]) is given by 

Hx =Y. 



+ /^o- (8) 
Similarly, the third example 

= aax, = iba-y, = ccr^ (9) 
belongs to the Ty phase and its parent Hamiltonian is 



+ (i - a^b^^)Sz,iSz,j + b'^c^-f)SyzaSyz,j 

+ (^ - a^C^'-f)SxzsSxz,j + + a^b'^l)SxyASxy,j 

+ /io- (10) 
The last example 

A'^ = oCTj;, = bay, A^ = icaz (11) 
belongs to the phase with its parent Hamiltonian given 

by 



E 



(^ - b'^c^l)Sx,iSx,j + ^^<^^l)^y,i^y,i 



,1 
.1 



'(^ + a%'^^)Sz,iSzj + (t + b'^c^'j)Syz,tSyz.j 



+ (- + a c 'y)Sxz,iSxz,] + (- - a 6 "f)Sxy,'i.Sxy,j 

+ ho. (12) 

Above we have given four special models that belong 
to different SPT phases. In the next subsection, we will 
show that if one keeps the D2h symmetry, phase transi- 
tions must happen when connecting these models. 



B. transitions between different SPT phases 

In order to justify the four SPT phase, we will use 
numerical method to study more general Hamiltonians. 
The method we adopt is one version of the tensor renor- 
malization group approach developed in ID by G. VidaP'^ 
and later generalized to 2D by T. Xiang et.al^ In this 
method, the ground state is approximated by a MPS. For 
an arbitrarily initialized state, we can act the (infinites- 
imal) imaginary time evolution operator U{St) — e"^"^"^ 
for infinite times, finally obtaining the fixed point matrix 




-E'(>1) 
-E"(tl) 
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FIG. 2: (Color online) (a) The energy curve E{ri), (b) The 
first and second derivatives of E{t)). All these curves are 
smooth, indicating that these states are in the same phase. 



A™. If the dimension D of A™ is not too small, the cor- 
responding MPS is very close to the true ground state. 
In our numerical calculation, we set D = 16. In ID, 
the ground state energy, correlation functions, density 
matrix, and entanglement spectrum can be calculated 
directly from the matrix A™. 

Noticing that /iq is a common term in the four exactly 
solvable models, which indicates that it is unimportant 
and can be dropped. This can be numerically verified. 
For this purpose, we add a perturbation to the models, 
such as 



H{ri) = Hx - rjho 



(13) 



where 77 S [0, 1]. As shown in Fig. [21 the ground state en- 
ergy E{ri) and its derivatives E'{r]), E"{ri) are all smooth 
functions, indicating that all the Hamiltonians H(r]) be- 
long to the same phase. Using the same method, one can 
also check that the Hamiltonians (with D2h symmetry) 
in the vicinity of an exactly solvable model fall in the 
same phase. For instance, the Heisenberg model and Ho 
in ^ are in the same phase. 

Now a question is whether the ground states of dif- 
ferent exactly solvable models can be smoothly trans- 
formed into each other. For this end, we consider a more 
realistic model ([T|) which connects two exactly solvable 
models, such as Hq and Hx- We are interested in the 
anti-ferromagnetic cases and will focus on the parame- 
ter region 9,4> & [0, §]• The point (^,0) is the Heisen- 
berg model. From the result of the last paragraph, the 
Heisenberg model is in the same phase as ([6]), and sim- 
ilarly (f , f ) is in the same phase as ([5]). If these two 
points cannot be smoothly connected (i.e., if gap closing 
or level crossing will unavoidably happen), then ^ and 
(|8]) belong to different phases. 

Using the tensor RG method, we can calculate the 
ground state energy of ^ and the phase diagram is 
shown in Fig. [T] When 9 is less then 0.2l7r, the ground 
state is Neel ordered. When 9 increases, a second order 
phase transition occurs and we enter the SPT phases. 
Fig. 13] shows the data of energy curve with fixed (j) = 
(and (f) — ^) and its first and second derivatives , which 
illustrate this transition. 

The region cj) G [0, |-) belongs to the Tq phase and 
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To and Tx phases the entanglement spectrum is doubly 
degenerate (see Fig. |3(b)| . This shows that the Tq and 
Tx phases are indeed nontrivial. Similar to the model 
the phase transition between Tq and Ty or between 
To and T^ can also be illustrated. 

Now we will show that first order phase transition also 
exists between any two of Tx,Ty,Tz. As an example, we 
consider the model that contains the transition between 
Tx and T^ phases: 

i 

~\~Sxy,iSxy,j) + Sai6(ySz,iSz^j + Syz^iSyzj)]- 

(15) 



When = tan i , the above Hamiltonian is in the same 
phase as Hx as shown in (fT3)) . and when 6 = tan^^ 5, it 
is in the same phase as Hz- The ground state energy of 
([T5| as a function of 9 can be obtained using the tensor 
RG method, and the result is shown in Fig. 01 A first 



(b) 



order transition at 9 
similar to 
Ei^~9). 



manifests itself. For the reason 
l4|) . the model also has a symmetry E{9) = 



FIG. 3: (Color online) (a) The energy curve E{9) with fixed 
<j!> = -I (and (f> — ^ ). (b) The line with circles is the first 
order derivative and the line with crosses is the second order 
derivative of E{6); the doted line is the difference tween the 
two biggest weights of the entanglement spectrum Ai — A2. 
Since E' (9) is continuous but E"{9) is not, the transition is 
second order. The degeneracy of Ai and A2 indicates that the 
To (and Tx) phase is nontrivial. 



(j> e (f , f ] belongs to the Tx phase. A first order phase 
transition between them happens at cj) — j. 

The first order phase transition occurs exactly at cj) — 
J . This is because the whole phase diagram is symmetric 
above and below the line (f> — j. This symmetry can be 
seen in the Hamiltonian. Notice that under a unitary 
matrix ?7, we get 



WSxU^Sx, WSyU^-Sxz, WSzU ^Sxy,{^4^) 



where U 

Hamiltonian 
H{9,^- 



This means that 



1 

,756-^1 ^e- 

© satisfies iU^(^U^)^H{9,^){ll^m^) = 
which yields E{9,(j)) = E{9,^ - 0), and 
their ground states are transformed by the (local) unitary 
transformation Yii ®Ui. But this unitary transformation 
is not invariant under time reversal T, so the behavior of 
the ground state also changes under time reversal. Re- 
sultantly, the state after the transformation belongs to a 
different phase. 

The entanglement spectrum of the ground states can 
also be obtained programmatically. We find that in both 




FIG. 4: (Color online) The ground state energy for model 
(|15|) . A first order phase transition is obvious. 

From the above analysis, we can conclude that the four 
exactly solvable models really stand for four distinct SPT 
phases. All these SPT phases are protected by the D2h 
symmetry. As will be shown in section IV, no more SPT 
phases exist for 5 = 1 models with D2h symmetry. Fur- 
thermore, the Eqs. ([T]) and (fTS]) show that these SPT 
phases can be obtained by much simpler Hamiltonians 
which is hopefully realized experimentally. 

Now an interesting question arises, how to distinguish 
these SPT phases in a practical way? It is impossible to 
distinguish these phases by linear response in the bulk 
since it is gapped. However, the end 'spins' localized at 
open boundaries may have different behaviors in different 
SPT phases. In the next section, we will propose an 
experimental method to detect each SPT phase. 
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III. DISTINGUISHING DIFFERENT SPT 
PHASES 



Xx(T),Xy{T), Xz{T) still holds, except that it is no longer 
isotropic. 



We expect to distinguish the four SPT phases through 
their different physical properties. Experimentally, all 
measurable physical quantities are response functions, or 
susceptibilities. So we need to add small perturbations 
and expect that (the 'end spins' of) different phases have 
different responses. The simplest perturbation for spin 
system is magnetic field H' = g^/isB-S, here S = Sj, 
(7i is the Lande factor and fig is the Bohr magneton. We 
will study the linear response to small B. 

Since the states in the same phase have the same 
universal properties, we will focus on the exactly solv- 
able models first. For simplicity, we consider the AKLT 
model, namely, Hq with a = b — c — 1. In the matrix 
product state picture, the physical S = 1 spin is divided 
into two J = 1/2 virtual spins. In the AKLT state, the 
virtual spins pair into singlets (called valence bonds) on 
each link between neighboring sites. Under open bound- 
ary condition, a free J = 1/2 spin at each end remains 
unpaired. The two end spins account for the exactly 
four-fold degeneracy of the ground states. In this pic- 
ture, it's easy to calculate the total spin in the ground 
state Hilbert space. The singlets in the bulk have no 
contributions to S, only the two end spins ji and j2 con- 
tribute and resultantly S = ji -I- j2- In this sense, the end 
spins can be considered as impurity spins of a paramag- 
netic material. Since the total spin of an open chain is 
^(815 = 001, we expect that the eigenvalue of Sx,Sy, Sz 
should be 1, —1, 0, 0. This can be verified by exact diag- 
onalizing a short chain. We denote these four degenerate 
ground states as 1-01), |'02), |'03)j l'04)- Then the matrix 
element of Sm in the ground state Hilbert space is given 
by 



s: 



a/3 



(V'al^mlV'/j), a,/3 = 1,2,3,4 



(16) 



The eigenvalues of the matrices (S"^) are exactly 
1, — 1, 0, and these values are independent of the length 
of the chain. Thus a small magnetic field along any di- 
rection H' = gLtJ-BBxSx or gLfJ-sBySy or ghtJ-sB^Sz will 
split the ground state degeneracy and give rise to a finite 
magnetization. 

At finite temperature, the susceptibility satisfies the 
Curie law and is given by21 



(17) 



where N is the number of end spins, g ~ \J J{J + 1)(7l, 
J = 1/2 and gL is the Lande g factor. If the spin-1 chains 
in the sample are broken into long separate segments, 
then N can be a considerable number. We also note 
that, in real samples, the susceptibility also contains a 
temperature independent part coming from the bulk. 

We see that, for the AKLT model, the spin suscepti- 
bility diverges at low temperature along all directions. 
For a general model in the Tq phase, the divergence of 




3 4 5 6 

Length of the chain (L) 

FIG. 5: (Color online) The eigenvalues of Sy and Sz are 
s, — s, 0, 0. The magnitude of s exponentially decays with the 
length of the chain in Tx phase. The dashed line is an expo- 
nential fit. The results are obtained by exact diagonalization 
and we only calculate up to seven sites. 

However, in phase T^, the end 'spins' have absolutely 
different physical properties. We consider the model 
in ([8]), and set a = 5 = c = 1. Then we calculate the 
eigenvalues of operators , Sy and in the ground state 
Hilbert space as before. We find that the eigenvalues of 
{S"^)x are still 1,-1,0,0, meaning that along x direc- 
tion the spin- 1/2 end spins still exist and Xx{T) diverges 
at T = 0. The eigenvalues of {S'^>^)y and {S'^'^)^ also 
have the structure s,—s, 0,0, but the magnitude of the 
nonzero eigenvalues s exponentially decay to zero with 
the increasing of the length of the chain (see Fig. [5]). This 
means that in y- and z-directions, there are no free spins 
coupled to the magnetic field. In appendix |D] we will 
show that this property is determined by the projective 
representation carried by the virtual 'spins'. In this case, 
Xy{T) and XziT) are given by ^ with 3^,5^ w 0. The 
result that Xx{T) follows Curie law and Xy{T), Xz{T) 
has effective gy,^^ « is a universal property of all the 
models in the phase. 

Similarly, one can check that only Xy {T) follows Curie 
law in Ty phase with the usual gy ^ J( J + l)gL 
[gxT gz ~ 0), and similarly only Xz{T) follows Curie law in 
Tz phase with the usual g^ « ^y J{J + l)gL {gx,gy ~ 0). 
Therefore, by measuring the temperature dependence of 
susceptibility and the effective g in x,y, z directions, we 
are able to distinguish the four SPT phases. 



IV. PROJECTIVE REPRESENTATIONS AND 
SPT PHASES 

In previous sections, we have given four SPT phases 
of the model and studied their physical properties. 
In this section, we will explain how the D2h symmetry 
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supports the existence of these phases. Then we will 
discuss other possible SPT phases of spin systems with 
D2h symmetry. 

The ground state of a gapped phase is written as ([3]) . If 
we require that the ground state MPS be invariant under 
the symmetry group D2h, namely, g\ip) = \ip) {g G D2h), 
then under the action of the symmetry group the matrix 
A™ must vary in the following way. 

1) g is unitary, g G {E, R^, Ry^R^,}, 

J2 <9)mm'A"^' = M{g)^A^Mig); (18) 

m' 

2) g is anti-unitary, g G {T, RxT^ RyT^ R^T}^ 

J2 u{gU^.KA^' = M{g)^A^M{g). (19) 

m' 

Here u{g) and M{g) are representations of the symmetry 
group D2h- The matrices u{g) satisfy the same multi- 
plication law of D2h group and is called a linear rep- 
resentation. The physical spin freedoms are linear rep- 
resentations of D2h- M{g) and M{g)e^^ are equivalent 
and belong to the same presentation. Up to a phase fac- 
tor (which depends on the group elements), M{g) satisfy 
the multiplication law of I?2/u and this kind of presen- 
tation is called a projective representation. The virtual 
'spins' (and the end states) are projective representations 
of D2h- More knowle dge a bout projective representation 
can be found in Ref. [aizil. 

The D2h group has eight 1-d linear representations (see 
Tab. IVI[) and eight (and only eight) classes of projective 
representations labeled by (111), (11-1), (1-11), (1-1-1), (- 
111), (-11-1), (-1-11) and (-1-1-1) (see Tab.|VlIl)- The first 
class of projective representation (111) is the eight 1-d 
linear representations, which are trivial (an example of 
the the corresponding trivial phase is the case 63 — > 00). 
The other seven projective representations are 2-d and 
nontrivial. Since states supporting different projective 
representations (or virtual 'spins') cannot be smoothly 
transformed into each other, each projective representa- 
tion corresponds to a SPT phase. This means that there 
should be at least seven different nontrivial SPT phases 
for spin systems respecting D2/1 symmetry. (In fact, when 
considering different a(g) there are more nontrivial SPT 
phases for spin systems respecting D2h symmetry.—) 

How can we obtain the projective representations? 
Mathematically, finding the projective representations of 
a group G is equivalent to find the linear representations 
of its cover group, which is a central extension of G and 
is called representation group R{G)^ The representa- 
tion group R{D2h) is available in literatures^ so we can 
calculate the matrix elements of all the projective repre- 
sentations of Z?2/i (see Tab. IVII|) . 

Once the matrices of the projective representations are 
obtained, we can calculate the CG coefficients for de- 
composing the direct product of two projective repre- 
sentations. From the CG coefficients, we can construct 
sMP states and their parent Hamiltonians^ The models 



(l6)).((8t. (|T0)) and are constructed accordingly and cor- 
respond to the (-1-1-1), (-1-11), (-11-1), (-111) representa- 
tions respectively. From these models we can know what 
kind of interactions are essential for each SPT phase. 

As shown in appendix [Bl the remaining three SPT 
phases of (1-11), (11-1), (1-1-1) cannot be realized for S — 
1 spin chains. The reason is that the physical freedom is 
not sufficient to support the direct product of two such 
projective representations. However, these phases might 
be realized in S* = 1 spin ladders or 5 = 2 models, and 
this will be our upcoming work. 

V. CONCLUSION AND DISCUSSION 

In summary, we have found four nontrivial SPT phases 
To,Tx,Ty,Tz of 5 = 1 spin chains which have on-site D2h 
symmetry. These SPT phases have similar properties 
as the usual Haldane phase, such as the bulk excitation 
gap, short-range correlations, existence of end 'spins', en- 
tanglement spectrum degeneracy. However, the different 
projective representations of the end spin under D2h in- 
dicate that they do belong to different phases. The SPT 
order that distinguishes them is the class of projective 
representations (or the group elements of the second co- 
homology iJ^ (1)2/1, C)) correspond to the ground states 
(or the matrices A™). 

We find that different SPT phases can be distinguished 
by experimental method. The magnetic susceptibilities 
Xx^XyiXz obey Curie law and diverge at zero tempera- 
ture. In To phase the effective (7-factors of the end spin 
have the usual values for magnetic field in x-, y-, and 
^-directions. But in Tx (or Ty or Tj) phase, the effective 
g sa J[J -\- l)gL has the usual value only for magnetic 
field in x-direction (or y-direction or z-direction). The 
effective g ~ [see eq. pT|) ] for magnetic field in the 
other two directions. 

From the seven nontrivial projective representations 
for D2/1 group, we constructed seven SPT phases. Four 
of them are discussed above. The other three may be 
realized in S* = 1 spin ladders or spin 5 = 2 models and 
are not discussed in the current paper. Some conclusion 
in this paper can be generalized to larger spin systems 
and higher dimensions. 
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Appendix A: Spin chian with D2 symmetry 

In appendix A, we will first study S — 1 spin systems 
with D2 symmetry. The same method can be applied to 
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D2h case. 



1. General Hamiltonian with D2 point group 
symmetry 

The point group D2 has only four elements, D2 = 
{E, Rx, Ry, Rz}- The multiplication table is shown in 
Tab.lll 



TABLE I: Multiplication table of D2 





E 


Rx 


Ry 


R. 


E 


E 


Rx 


Ry 


Rz 




Rx 


E 


R. 


Ry 


Ry 


Ry 


R. 


E 


Rx 


R. 


R. 


Ry 


Rx 


E 



It has four 1-d linear representations, whose matrix el- 
ements and bases of representations are shown in Tab.llll 
From quantum mechanics, we know that the 25' + 1 bases 
of integer spin-S' span an irreducible linear representation 
space of SO{3) group. This Hilbert space is reduced into 
a direct sum of 25* + 1 1-d irreducible linear representa- 
tion spaces of D2- For example, when 5 = 1 (a vector), 
the bases 

k> = ^(|-l>-|l», 
|y) = i=z(|-l> + |l)), 
= |0) 

form the B3, B2,Bi representations of D2 respectively. 



TABLE IL Linear representations of D2 





E 


Rx 


Ry 


R. 


bases or operators 


A 


1 


1 


1 


1 


10,0) 




q'Z q'Z q'2 
•> •> ^ z 


Bi 


1 


-1 


-1 


1 






Sxy 


B2 


1 


-1 


1 


-1 


|i,y> 


Sy 


Sxz 


B3 


1 


1 


-1 


-1 






Syz 



Here we focus on the S = 1 model with nearest neigh- 
bor interaction. The general Hamiltonian with D2 sym- 
metry is given by 

'^f2{Sy^iSxz,j "I" Sxz,iSyj) 

~^f3{Sz,iSxy.j + Sxy^iSzj)- (^1) 

where /i,/2,/3 are constants and Ho^^ is given in ©. 
The above Hamiltonian Hu^ also has translational sym- 
metry and spacial inversion symmetry. The additional 
/i , /2 , /s terms are odd under time reversal and break 
the T symmetry of Hjj^. To study the SPT phases, we 
need to obtain the projective representations of D2. 



2. Projective representation and CG coefficients of 

D2 

From Ref. HJ, determining the projective representa- 
tion of a point group G is equivalent to determining 
the linear representation of its representation group (s) 
R{G) (which cover G integer times). There are two 
non-isomorphism representation groups of D2, namely, 
Ri{D2) and i?2(^2), both of which have two generators 
P, Q and eight group elements. Their multiplication ta- 
bles are listed in Tabs. IHII and IIVI In the following, we 
will mainly discuss the covering group i?i(D2), and leave 
the discussion about i?2(£'2) to the end of this section. 



TABLE III: Multiplication table of Ri {D2). Notice that P* = 
Q* = E, P'^ ^ g2 and QP = P^Q. 





E 


p2 


p3 


P 


Q 


P'Q 


PQ 


P'^Q 


E 


E 


p'2 


p3 


P 


Q 


P'^Q 


PQ 


P^Q 


p2 


p2 


E 


p 


p3 


P'Q 


Q 


P^'Q 


PQ 


pi 


p3 


P 


p2 


E 


P^'Q 


PQ 


Q 


P^Q 


P 


P 


p3 


E 


p2 


PQ 


P^Q 


P^Q 


Q 


Q 


Q 


P^Q 


PQ 


P'^Q 


p2 


E 


P 


p3 


P^Q 


P^Q 


Q 


P^'Q 


PQ 


E 


p2 


p3 


P 


PQ 


PQ 


P^Q 


P^Q 


Q 


p3 


P 


p2 


E 


P^Q 


P^Q 


PQ 


Q 


P'Q 


P 


p3 


E 


p2 



TABLE IV: Multiplication table of ^2(1)2). Notice that P* = 
Q'^ = E and QP = P^Q. 





E 


p2 


p3 


P 


Q 


p'g 


PQ 


P'^Q 


E 


E 


p2 


p3 


P 


Q 


P'^Q 


PQ 


P^Q 


p2 


p2 


E 


P 


p3 


P^Q 


Q 


P'^Q 


PQ 


p3 


p3 


P 


p2 


E 


P^Q 


PQ 


Q 


P^Q 


P 


P 


p3 


E 


p2 


PQ 


P'^Q 


P^Q 


Q 


Q 


Q 


P'Q 


PQ 


P^Q 


E 


p2 


p3 


P 


P'^Q 


P^Q 


Q 


P^'Q 


PQ 


p2 


E 


P 


p3 


PQ 


PQ 


P^'Q 


P^Q 


Q 


P 


p3 


E 


p2 


P'^Q 


P^Q 


PQ 


Q 


P^Q 


p3 


P 


p2 


E 



To obtain all the irreducible representations, we only 
need to block diagonalize the canonical representation 
matrices of the two generators P and Q. In the canonical 
representation, the group space itself is also the represen- 
tation space. Each group element 171 is considered as an 
operator gi: 

91(92) =gi92- (A2) 

here 32 and 51 (?2 are two vectors in the representation 
space and gi becomes a matrix. 

The canonical representation matrices of the genera- 
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tors of i?i can be read from table IIIII 



P = 



/O 1 0\ 




/O 


1 





0\ 


1 




1 











1 













1 


1 










1 





1 


1 











1 




10 











1 




10 











\0 1 0/ 




\0 1 








0/ 



To simultaneously block diagonalize the above matrices, 
we need to identify the base vectors (or wave function) of 
each irreducible representation (these base vectors form a 
unitary matrix which block diagonalize P and Q simulta- 
neously). In quantum mechanics, we use good quantum 
numbers (eigenvalues of commuting quantities) to label 
different states. For example, \S, m) symbols a spin state, 
where S{S + 1) is the eigenvalue of the Casimir operator 
of SO{2>) group and m the eigenvalue of the Casimir op- 
erator of its subgroup 5'0(2). Similar method has been 
applied to the representation theory of groups. What 
we need to do is to find all the commuting quantities, or 
the complete set of commuting operators (CSCO)<^ 

The Casimir operators of discrete groups are their class 
operators. For i?i(I?2), there are five classes (hence there 
are five different irreducible linear representations), and 
the corresponding five class operators are given as below: 

C^{E, P\ P + P^, Q + P^Q, PQ + P^Q}. 

The class operators commute with each other and all 
other group elements. This set of class operators C is 
called CSCO-I in Ref. [2^ The eigenvalues of the class 
operators are greatly degenerate, which can only be used 
to distinguish different irreducible representations(IRs). 
To distinguish the bases in each IR, we can use the class 
operators of its subgroup (s). Group i?i has a cyclic sub- 
group 

C{s)^{E,P,P\P'}, 

each element forms a class. The set of class operators 
of the subgroup is written as C(s). The operator-set 
(C, C(s)) is called CSCO-II, which can be used to dis- 
tinguish all the bases if every IR occurs only once in the 
reduced canonical representation. 

However, in the reduced canonical representation, a d- 
dimensional representation occurs d times and they have 
the same eigenvalues for CSCO-II. To lift this degeneracy, 
we need more commuting operators. Fortunately, we can 
use the class operators of the 'intrinsic group' Ri, whose 
group elements are defined as follows 



31(32) = 32fll- 



(A3) 



Notice that g commutes with g defined in (IA2p . The class 
operators of Ri are identical to those of C = C. The 
set of class operators for the intrinsic subgroup 



is noted as C{s). The eigenvalues of C{s) provide a dif- 
ferent set of 'quantum numbers' to each identical IR. 

Now we obtain the complete set of class operators 
{C,Cis),C{s)), which is called CSCO-III. The common 
eigenvectors of the operators in CSCO-III are the or- 
thonormal bases of the irreducible representations, and 
each eigenvector has a unique 'quantum number'. 

To obtain the bases, we need to simultaneously diago- 
nalize all the operators in CSCO-III and get their eigen- 
vectors. Actually, we only need a few of these operators, 
for example, we can choose Q + P^Q in C, P in C(s) and 
P in C{s). The matrices of these operators of Ri{D2) 
are given below: 



Q 



p = 




















1 


1 

























1 


1 
































1 


1 


























1 


1 


P Q = 




1 


1 


























1 


1 
































1 


1 




















Vo 





1 


1 











oy 







1 














^\ 

















1 























1 


























1 


















































1 


1 
























1 























1 

















Vo 














1 





0^ 













1 














^\ 

















1 























1 


























1 















































1 





; 



























1 























1 














Vo 











1 








oy 









Practically, we can diagonalize a linear combination 
O = {Q + P^Q) + o,P + bP, where a, h are arbitrary 
constants ensuring that all the eigenvalues of O are non- 
degenerate. From the non-degenerate eigenvectors (col- 
umn vectors) of O, we obtain an unitary matrix U : 



U = 





1 


1 


1 


































V8 


\/8 


V8 



1 
2 










1 
2 








1 

2 


-¥ 








¥ 


¥ 








-¥ 







-¥ 








¥ 


¥ 








1 

2 


1 
2 








1 
2 


1 
2 


) 



(A4) 



c(s) = {£;,p,p2^p3} 



The matrix U is the transformation that block diagonal- 
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izes P and Q simultaneously: 



(the column vectors are just the CG coefficients): 



U^QU 



/I 














^\ 





1 














n 


n 


_1 


n 


n n 

u u u 


n 


n 
U 


u 


u 


1 

— 1 


n n n 
u u u 


n 
u 


A 

U 


u 


u 


u 


— I u u 


u 


u 


u 


n 
u 


u 


U i u 


u 


u 


u 


n 
u 


u 


u u — t 


A 

u 


\0 


u 


u 


u 


n n n 
u u u 


V 


/I 











0\ 







-1 



















1 



















-1 



















i 
















i 
















i 




Vo 











i 0^ 





(A5) 



(A6) 



There are four 1-d IRs and one 2-d IR (which occurs 
twice) in the reduced canonical representation: 



-i 
i 



i 

1 



(A7) 



So there are totally five independent representations. 

The eight elements in Ri can be projected onto D2, as 
shown in Tab. [V] And the linear representations of Ri 
correspond to the projective representations of D2. The 
four 1-d IRs correspond to the linear IRs of D2, and the 
2-d IR stands for a nontrivial projective IR of D2. Up to 
a phase factor, these 2-d matrices are the 180° rotation 
operators of a spin with J = 1/2 (which is a projective 
IR of SO (3) group). 



TABLE V: Projection from Ri{D2) to D2 



Ri{D2) 


E 

p2 


P 

p3 


Q 
P'^Q 


PQ 
P^Q 


D2 


E 


R. 


Rx 


Ry 


rotation tt of J = 1/2 










(up to a phase factor) 


I 






iay 



Now let's look at the direct product of the projective 
IRs of Z?2. For the 1-d linear IRs, the direct product are 
still 1-d IRs, which satisfy the following law: 

Ax Bi^ Bi, Ax B2 = B2, Ax B3= B3, 
BixB2 = S3, BixB3^ B2, B2 X S3 = Bi. 

The direct product of 1-d and 2-d IRs are still 2-d projec- 
tive IRs of D2. The direct product of two 2-d projective 
IRs is interesting. It reduces to four 1-d linear IRs. Us- 
ing the CSCO-II, we can diagonalize the 4x4 matrices 
P2 (8) P2 and Q2 <E) Q2 with the following unitary matrix 



U4 = 



V2 



I 

V2 










1 


1 




\/2 





\ 




V 


V 


i 





/ 







0\ 




-1 


0] 


(• 





1 







1/ 







0\ 




-1 





(= 







Vo 





1/ 



Ul{P2 <E) P2)U4 = 



UI{Q2®Q2)Ua 



These CG coefficients of the projective IRs of D2 are 
analogous to the decoupling of the direct product of two 
spins with J=l/2, — \ ®Q except the 3-d IR of 

spin-1 becomes a direct sum of three 1-d IRS. If we label 
the bases of the 2-d projective representation of D2 as 
I t), I I), then the CG coefficients are given as: 



\X) = i=(| UU) - I tlt2)), 

|y) - -^(1 Ills) + I tit2)), 
|z) = -i=(| tiiz) + I iits)), 
|singlet) = -i=(|ti;2)-Uit2)), (A8) 



Repeating the above procedure, we obtain the IRs of 
R2{D2)- The four 1-d IRs are the same as that of _Ri(Z?2), 
while the 2-d IR is given as: 



P' - 
^2 — 



-i 



i^2 ~ \ \ Q 



(A9) 



The above representation and Eq. (jA7l) differ only by a 
gauge transformation = P2 and Q'2 = iQ2, so they 
belong to the same projective representation of D2. The 
CG coefficients for the 2-d IRs are obtained easily: 



\x) 


= ^(1 ;i;2) + 


1 tlt2)), 


\y) 


= ^(1 ;i;2) - 


Itlt2)), 




^^(It42)- 


1 ;it2)), 



Isinglet) = ^(1 tiU) + I ;it2)). 



(AlO) 
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3. sMP state with D2 symmetry and its parent 
Hamiltonian 

Before studying the model with D2 symmetry, let's 
review the 5 = 1 AKLT model^^ (which has S'0(3) sym- 
metry) first. The AKLT state is a sMP state given by 
= ax,Ay = (Ty,A^ = az- The A™ matrices are two 
by two, meaning that the physical spin S" = 1 is viewed 
as symmetric combination of two J — \/2 virtual spins 
(essentially projective representations of 50(3)). Alter- 
natively, we can write the state as 

=Tt{WiW2...Wn), (All) 

where Wi = A^\x)i + Ay\y)i + A''\z)i. According to 
Ref. [2^ the matrices A™ of a sMP state can be obtained 
by 

^™ = B'^(C")*, (A12) 

where B is the CG coefficient combining two virtual spins 
into a singlet |0, 0) — i?mim2 I^j \^ "^2), and C" is the 
CG coefficient combining two virtual spins into a triplet 

|l,w) = C^^^Ji, mi; i, ma). 

Now we can generalize this formalism to the D2 case, 
where the three states of S* = 1 become a direct sum 
of three IRs of D2. D2 group has a 2-dimensional non- 
trivial projective representation, and the direct product 
of two such projective IRs can be reduced using the CG 
coefficients {B and C'^''^) in Eq. (|l8|. Similar to the 
50(3) case, we can consider the two 2-d projective IRs 
as 'virtual spins'. From eq. (IA12[) . we can construct the 
following matrix (similar sMP state has been studied in 
Ref.|27D 

W = aax\x) +bay\y) + caz\z), (A13) 

where a, &, c are arbitrary nonzero complex constants. 
The corresponding sMP state is given by |<^) = 
Tv{WiW2---Wn), which is invariant under the group D2. 
Notice that the CG coefficients in Eq. (jAlOp give the 
same sMP state (up to some gauge transformations). No- 
tice also that (|A13|) is different from ©, 0, (l9|) or ([TT1). 
If a, 6 or c are to be arbitrary complex numbers, it is not 
invariant under T . 

The above sMP state is injective, and the parent 
Hamiltonian can be obtained by projection operators. 



We consider a block containing two spins, the four matrix 
elements of WiWi+i span a 4-dimensional Hilbert space. 
Suppose the orthonormal bases are |'(/'i,2,3,4)i, then we 
can construct a projector 

4 

P, = l-Y,\^c.){il>cU (A14) 

a = l 

and the Hamiltonian H = J^i Pi- It can be easily checked 
that the sMP state is the unique ground state of this 
Hamiltonian. 

The projector Pi is a nine by nine matrix that can 
be written in forms of spin operators. Notice that any 
Hermitian operator of site i,j can be expanded by the 81 
generators of L''(9) = U{3)i (g) U{3)j, i.e., Xai^pj {01, (3 — 
1, 9). So, we have 

9 

^'i = ^ £.al3>^ai^Pj, (A15) 
a,/3 = l 

where are constants. Further, the generators of U{3) 
can be written as polynomials of spin operators. 

Ai = {Sx + Sxz)/V2, 
Ae = {Sx - Sxz)/V2, 

A2 = (Sy + Syz)/V2, 

Xj = {Sy - Syz)/V2, 

A4 »5*^ Sy , 

A5 — Sxy-, 

A3 = (5,+352)/2-/, 

As = {?,Sz - 352 + 2/)/2V3, 

where 5„i„ = 5,„5„-f 5„5„i, {m,n = x,y,z) and Ai As 
are the Gellmann matrices of SU (3) generators. Finally, 
we can write the Hamiltonian in forms of spin operators. 
For simplicity, we first assume a, &, c are real numbers, 
then the Hamiltonian is given in which is invariant 
under T. The T symmetry goes away when a, 6 or c 
becomes an arbitrary complex number. For instance, if 
a — > ae*^, then the Hamiltonian ^ becomes 



H = 



E 



(7 + 



,1 cos29a^c^/2 



,4 )Sx,^Sx,J + ( 4 + ^4 + 54 + ^4 )Sv,rSy^J ^4 



sin26la2cV2 



fe4 



1 <MsWa^lPj2 sin 26*02672 

+ ^7 + „4 1^4 ^ ^4 ^-^^^^-^^J + " 



&4 



l__6V/2_ 1 cos2ea\^/2 

^4j_54 , „i)^V^,^^V^,3 + \ A „4 , ;,4 , ^ )'^xz,ii:}xz,3 



{Sz,iSxyJ + Sxy,iSzj) + { ^ ^4 

ho 



1 cos26la2fe72 



64 



64 



(A16) 
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When sin 20 7^ above Hamiltonian does not have T 
symmetry. 

Varying the values of a, h, c, we can transform the 
ground state of the above Hamihonian into that of the 
AKLT model smoothly without breaking D2 symmetry. 
This means that above sMP state also belongs to the Hal- 
dane phase. In appendix [B] we will consider the models 
with additional time reversal symmetry. 

Appendix B: Spin Chain with D2h symmetry 

In the last section we have studied the spin chain with 
on-site D2 symmetry. Now we consider a 5* = 1 spin 
chain with additional spin-inversion (or time-reversal) 
symmetry. The complete on-site symmetry now becomes 
D2h = {E,R,,Ry,R„T,R^T,RyT,R,T}. It has eight 
1-d linear real IRs, as listed in Tab. IVII Notice the time 
reversal operator T = e~'^^^'"K is anti- unitary, so the 
states |m) and i\m) (m = x,y,z) belong to different lin- 
ear representations, the former is odd under T and is 
noted by index u, and the latter is even under T as noted 
by g. So we need to introduce six bases \y), \z) and 
i\x),i\y),i\z). To construct a sMP state, at least one of 
the pair \x), i\x) (and also the pairs \y), i\y) and \z), i\z)) 
should be present in the physical bases. 



Here all the coefhcients are chosen to be real. 

Now we construct sMP states from the CG coefficients 
Eqs. dHT]), (lB2|) and ((M2)) . Since aU the CG coefficients 



TABLE VI: Linear representations of D2h 





E 


Rx 


n 

IXy 


Rz 


T 


RxT 




RzT 


bases 


Operators 


A„ 


1 


1 


1 


1 


1 


1 


1 


1 


10 0) 


5^ S'^ S'^ 

^x') ^yi ^ z 


Big 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


i\l,z) 


Sxy 


B2g 


1 


-1 


1 


-1 


1 


-1 


1 


-1 




Sxz 


B'ig 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


i 1, x) 


Syz 


Au 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


i\0,0) 




Blu 


1 


-1 


-1 


1 


-1 


1 


1 


-1 




Sz 


B2u 


1 


-1 


1 


-1 


-1 


1 


-1 


1 


|i,y> 


Sy 


Bzu 


1 


1 


-1 


-1 


-1 


-1 


1 


1 




Sx 



To obtain the projective IRs of -D2/1, we need to 
study the linear IRs of the representation group R{D2h), 
which also has three generators P,Q,R (corresponding 
to R,,R^,T) satisfying P'^ = Q'^ ^ R'^ = E and 
P^Q = QP, Q^R = RQ, R^P = PRJ^ The total number 
of elements in R(D2h) is 64. It has 8 1-d representations 
(corresponding to the 8 linear IRs of D2h) and 14 2-d 
representations (corresponding to the 7 classes of projec- 
tive IRs of D2h)- To obtain the IRs of R{D2h), we only 
need to know the representation matrix of the three gen- 
erators P,Q,R. Using the same method given in the last 
section, we obtain all the IRs of R{D2h) (see Tab. IVIip . 

Now we give the CG coefficients that reduce the direct 
product of two projective IRs into direct sum of linear 
IRs of D2h- 



(Bl) 



(B2) 



are real, the constructed matrices A™ = B'^{C"^)* are 
also real (here B = C'^", m = Big, Biu-.-B^u), and are 
invariant under the anti-unitary operator K. However, 



I 

Ei^Ei=E2(g>E2=Ag(S Big © A„ © = a^, C^'" = a„ C^" ^ ioy, C^^" = /; 

E-i®E-i=E^®E^ = Ag® B-ig © A„ © Bg^; C^^ = /, C^^^ = ioy, C^- ^ g„ C^'" ^ a^- 
E5<»E5=Ee<»Ee^Ag® Big © B2g © ^33; = 'Jx, C^'' - i<Jy, C^'' = <Jz, C^'" = I; 

Et (E) Er = Es (E) Es ^ Big © Bsg © -Bi„ © Bg^; C^^' = a„ C^^' = tay, C^'- = I, C^'- = a,; 
Eg(g)EQ= Eio (E> Eio = Big © Bag © A„ © ^3^; C^" = a,, C^'- = I, C^'^ = i<jy, C^^' = a,; 
Ell «) £^11 = E12 (g> E12 = B2g © Bsg © Au © Blu; C"^" = icTy, C^^" ^ cr^, C^^' = /, C^^« = a,; 
Ei3 <E) Ei3 = Eli ® Eu ^Ag® B2g © Bi„ © B^u] C^" = i(Ty: = /, C^^" = a„ C^^" = a,; 

Ei®E2 = B2g © B3g © B2u ® Bsu, C^'' = <Jx, C^'^ = a,, C^'" = tag, C^'" = J; 
Es(g)E4^ Big © B2g © Blu ® B2u; C^'' = /, C^'' = i<Jy, C^'" = <Jz, C^'" = <Jx; 
E5(g,Ee^ Au © Blu ® B2u ® Bs^, C^" = <Jx, C^'" = t<Jy, C^^" = a,, C^'- = /; 
Er(g)Es^Ag® B2g © A„ © B2«; C^" = cr,, C^^' = iay, C^" = /, C^^" = cr^; 
Eg Eio = Ag® Bsg © Blu © B2u; C^' = <7x, C^'' = /, C^^" = i(Ty, C^^" = a,; 
Ell <E) E12 = Ag® Big © B2u © S3„; C^^ = idy, C^i» = CT„ C^^" = /, C^^" = a,; 
Ei3 Eu - Big © Bsg © A„ © B2u; C^" = i<yy, C^^" = /, C^^' = a,, C^^' = a,. 

I 
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TABLE VII: Projective representations of D2h- The numbers q,/3,7 are obtained hy a — , j3 — , 7 = i?^. The three 
generators P,Q,R of R(D2h) will project to Rz,Rx,T of D2h, respectively. 







Q{Rx) 


R(T) 


a = P^, 3 = 0'^, 7 = R^ 


A 




1 


1 


1 




D 

-D19 




1 


-1 


1 




02g 




-1 


-1 


1 




D 

mg 




-1 


1 


1 


111 


A 




1 


1 


-1 








1 


-1 


-1 




Rn 




-1 


-1 


-1 








-1 


1 


-1 








I 


iaz 


ay 


1-11 




R.> 


-I 


iaz 


ay 




^3 






I 


iay 


11-1 


rp — p 


Ri 

-Dig 




-I 


iay 










ax 


I 


-111 






iaz 


ax 


-I 




Et 






iaz 


ia^ 


1 -1 -1 


E8^E7(g 


B19 


(Tz 


-iaz 


iax 




Eg 




iaz 


ax 


iax 


-1 1 -1 


Eio = Eg ( 




icTz 


ax 


-iax 




Ell 




iuz 


iax 


az 


-1 -1 1 


E12 ~ Ell c 


3 -Bag 


ia-z 


iax 


-az 




Eis 




irjz 


iax 


iay 


-1 -1 -1 


Ei4 — Eig 


® ^„ 


iaz 


iax 


-iay 





the bases \Big) — i\z), \B2g) — i\y) or l-Bsg) = i\z) con- 
tain a factor i, this factor i may be combined with A"^ 
when writing the matrix W = X^m^™!'^)- 
nition of depends on the choice of base. If we choose 
|m) = \x), \y), \z) as the physical bases, then A"^ will 
absorb the factor i (if existent) and may be either real 
or purely imaginary. This convention is adopted in the 
main part of this paper. On the other hand, if we just 
choose |m) = j-Bi^), |-Bi„), ...|-B3„) as the physical bases 
(and forget about the factor that some bases, such as Big 
and -Bim, are linearly dependent), then all the matrices 
A™ are real. In the following discussion, we will adopt 
the second convention. 

Notice that the combinations -B5 -E5, Eg (E) Eiq, 
Ell (8) E12, Ei3 (g) Ei3 contain all the bases of 5 = 1 
(l-Bi), I-B2), I -83)) and the singlet state (|^g)), we can con- 
struct sMP state using these combinations. We will study 
them case by case. 

1) E5 (8) E5 

Up to an overall phase, the local matrix W is given by 
W = a<7x\x) -\- ibay\y) -t- c<7z\z), here a, b, c are real num- 
bers. The Hamiltonian can be constructed using the 
method given in appendix |^ and the result is given in 

m- 

2) Eg Eio 

Up to an overall phase, the local matrix W is given by 
W = aax\x) -f b(7y\y) + iccr^jz), and the Hamiltonian is 
shown in (fT^ . 

3) £;ii (8) E12 

Up to an overall phase, the local matrix W is given by 
W = iaax\x) -\- bay\y) -\- caz\z), and the Hamiltonian is 



given in 

4)Ei3 (g) Ei3 

The local matrix W is given by = aax\x) -f bay\y) -\- 
caz\z), and the Hamiltonian is given in ([5]). 

With the -D2/1 symmetry kept, the ground states of the 
above four exactly solvable models cannot be smoothly 
transformed into each other, which indicates they belong 
to different SPT phases (see section HI). 

According to Ref. H, there should be seven SPT phases 
since there are seven classes of projective representations. 

However, in the other three projective IRs, the reduced 
Hilbert space of the direct product of two virtual 'spins' 
only contains one of the three bases for the physical S = 1 
states (notice that the singlet \Ag) is necessary to con- 
struct a sMP state), which means that these three SPT 
phases cannot be realized in S" = 1 systems. 

Appendix C: Invariance of the sMP state under 
symmetry group 

Firstly, we assume that all the operators of the sym- 
metry group G are unitary. The CG-coefficients (of the 
representation group) are defined as 

a,l3 

Isinglet) = EB„^|a,/3) (CI) 

a,l3 

where |m) belong to nontrivial linear IRs and | singlet) is 
a trivial linear IR, a, /3 are bases of some 2-d projective 
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IR. We will show that the sMP state given by (jA12[) is in- 
variant under the representation group i?(G) (and hence 
the symmetry group G). Suppose that is a group ele- 
ment oi R{G), and u{g)/N{g), M{g) is the representation 
matrix for the physical spin/virtual 'spins', then 



I to) = u„i'„i\m ), 
\a) = A^a'ala'), 

|/3) = M0,^|/3'), 



From Eqs. (|CTj) and ([02]), we obtain 



(C2) 



(C3) 



The complex conjugate of above equation is 

E"™™'(C^"')*=^*(C^")*Mt. (C4) 

m' 

Since the representation matrix u{g) (N{g),M(g)) 
is unitary, so the representation matrix of g^^ is 
[uig)]H[N{g)]^,[M{g)Y). Replacing g by g-^ in 
Eqs. (IC2|-(ICl|, we obtain 

Y.Umm'iC"'')* - Af^(C")*M. (C5) 

Similar to (|C3| . we also have 

B = NBM^, 
or equivalently — MB'^N^ . Thus we have 

M^B^ = B'^N^. (C6) 
From jCU and jCl]), we have 



\ m / 771, m' 

= J2 i?^"m'm(C")*|m') 
m,m' 

= Em^B'^(C"")*M|to) 

m 

= EmU"M|to). (C7) 



Above equation is nothing but (|4]), which indicates that 
the sMP state constructed by A™ = B^(C™)* is really 
invariant under the group R{G) (or equivalently, the sym- 
metry group G). 

Now we consider the case that some group elements, 
such as the time reversal operator T, of G are anti- 
unitary. Suppose that by properly choosing the phases 
of \a) and |/3), all the CG-coefficients B and C™ are real. 



In this case, the anti-unitary operators behave as uni- 
tary operators when acting on A^, and (jC7|) also holds 
for anti-unitary operators. 

To obtain the complete representation of the anti- 
unitary operators, we introduce an unitary transforma- 
tion to the bases \a), |/3) of the virtual 'spins' so that A™ 
transforms into complex matrix: 



\a) = Va'a\a') 



(C8) 



then 



= iVC"^U^)^,p,\a'P')=C':;^f,,\a'/3'), (C9) 

which gives C = VC™U^. Similarly, we have B' = 
VBU^. Since A"" = B^(C")*, so we get A" = 
UA™W. When an unitary operator g„ acts on A "^\m), 
(|C7P holds as expected: 

5„(Ea'"») = Ugu{J2A^\m})U^ 

m m 
rn 

= ^(M')U'"Af'|TO>, 

rn 

where M' = C/A/C/t. Now let us see what happens if an 
anti-unitary operator ga acts on A ™|to): 



5a E^'">^) = ff« UiJ2A"^\m))U^ 

\ m / \_ m 

= ^C/*MU™MC/^|to) 

rn 

= Y{M')^A"'M'\m), 



where M' = UMU^ = U{MK)W. Here we have used 
the factor that A"^ are real matrices. So (|C7p still holds 
for anti-unitary operators, except that the representation 
matrix M[ga) transforms into M'{ga) instead of M'{ga), 
or equivalently, M{ga) is replaced by M{ga)K . Thus for 
an anti-unitary operator g, we have 



L(g)K{A'^) = KAPA'^MK. 



(CIO) 



This result will be used in appendix [Pl 

Notice that to obtain a sMP state that is invariant 
under a symmetry group G containing anti-unitary op- 
erators, the only condition we require is that the CG co- 
efRcients B and C™ (for the unitary projective IRs of G) 
can be transformed into real numbers by choosing proper 
phases. 
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Appendix D: effective operators in the ground state 
Hilbert space 

From the projective representation, we can study the 
effective operator of a usual operator (which acts on the 
physical spin Hilbert space) on the ground state Hilbert 
space, or equivalently, the end 'spins'. Naturally, the 
usual operator and its effective operator should vary in 
the same way, or respect the same linear representation, 
under the group I?2h- So we will study the effective op- 
erators from the symmetry point of view. 

If the spin chain is long enough, the two end 'spins' are 
free (i.e., the interaction between them are neglectable). 
So we expect that the effective operators on the end 
'spins' are single-body operators instead of two-body in- 
teractions. Notice that the all the nontrivial projective 
representations of D2h are 2 dimensional, we have only 
three choices of the effective operators, the pauli matri- 
ces. We will study them one by one. 

Firstly, we study the Tq phase, which correspond to 
the projective IR (-1-1-1). Under symmetry operation g, 
the operator Om varies in the following way 

M{g)^d^M{g) - 77(5) wO™', (Dl) 

where M{g) is the projective IR for the end 'spin'. From 
the conclusion in appendix [C] and Tab. I VIII we get 
M{R,) = ia,, M{R^) = la,, M{T) = lUyK. rjig) is 
a linear representation of I?2h, which equals either 1 or 
— 1. Actually, ri{g) is the parity of Om under g. For in- 
stance, 77(T) = — 1 means that Om has odd parity under 
time reversal transformation and vice versa. After sim- 
ple algebra, we obtain the correspondence in table IVIIII 
the operators in the same column transform in the same 
way. 

TABLE VIII: Correspondence between physical operators and 
effective operators in To phase according to their transforma- 
tion property (parities) under D2h- 



linear IR ri{g) 




B'zu 


Blu 


operators Om 




Oy 




physical operators 




Sy 


Sz 



From above table, we find that am and Sm (w — 
X, y, z) have the same symmetry (or the same parity un- 
der symmetry operations), so the former can be consid- 
ered as the effective operator of the latter. Since am is the 
spin operator of the end spins, the system will response 
to weak external magnetic field (along any direction) ef- 
fectively through the end spins. 

However, things are different in Tx phase, which cor- 
responds to the projective IR (-1-11). From Tab. IVIH 
we can substitute M{Rz) = ia^, M{Rx) = ia^, M{T) = 
azK into (IDip and obtain the results in table lIXl Notice 
that the end 'spin' operator ay{az) do not have the same 
symmetry with that of Sy{Sz) because they have differ- 
ent time reversal parities. Since there are no single-body 



TABLE IX: Correspondence between physical operators and 
effective operators in Tx phase, according to their transfor- 
mation property (parities) under D2h- 



linear IR ri{g) 


B^u 


B2g 


Big 


operators Om 


Ox 


ay 




physical operators 


Sx 


Sxz 


Sxy 



effective operators correspond to Sy and Sz, the mod- 
els in the phase will not response to weak external 
magnetic fields along y- and z- directions. 

Similar results can be obtained for Ty and Tz phases 
and will not be repeated here. 
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